Negative-energy perturbations 
in cylindrical equilibria 
with a radial electric field 



G. N. Throumoulopoulos^ and D. Pfirsch 
Max-Planck-Institut fiir Plasmaphysik, EURATOM Association 

D-85748 Garching, Germany 

July 1997 



To be published in Physical Review E 



^Permanent affiliation: Section of Theoretical Physics, Physics Department, University of 
loannina GR 45 f 10 loannina, Greece 



Abstract 



The impact of an equilibrium radial electric field E on negative-energy per- 
turbations (NEPs) in cylindrical equilibria of magnetically confined plasmas is 
investigated within the framework of Maxwell-drift kinetic theory. It turns out 
that for wave vectors with a non-vanishing component parallel to the magnetic 
field the conditions for the existence of NEPs in equilibria with E = [G. N. 
Throumoulopoulos and D. Pfirsch, Phys. Rev. E 53, 2767 (1996)] remain valid, 
while the condition for the existence of perpendicular NEPs, which are found 
to be the most important perturbations, is modified. For |ej0| ~ Tj, a scaling 
which is satisfied in the edge region of magnetic confinement systems {(j) is the 
electrostatic potential), the impact of E on perpendicular NEPs depends on the 
value of Ti/Te, i.e., a) for T^/Te < /3c ~ P/{B^/8n) (P is the total plasma pres- 
sure) the electric field does not have any effect and b) for Ti/Te > /?c, a case 
which is of operational interest in magnetic confinement systems, the existence 
of perpendicular NEPs depends on e,^E, where e^, is the charge of the particle 
species u. In the latter case, for tokamaklike equilibria and H mode parameters 
pertaining to the plasma edge two regimes of NEPs exist. In the one of them the 
critical value 2/3 of 77^ = ^InTj/^ln A^j plays a role in the existence of ion NEPs, 
as in equilibria with E = 0, while a critical value of rje does not occur for the 
existence of electron NEPs. However, E has a "stabilizing" effect on both particle 
species in that (a) the portion of particles associated with NEPs (active parti- 
cles) is nearly independent of the plasma magnetic properties, i.e., it is nearly the 
same in a diamagnetic plasma and in a paramagnetic plasma, while in equilibria 
with E = this portion is much larger in a paramagnetic plasma than in a dia- 
magnetic plasma, and (b) the fraction of active particles can decrease from the 
plasma interior to the edge, e.g., for the case of electron NEPs in an equilibrium 
of a diamagnetic plasma, contrary to equilibria with E = 0. In particular, the 
fraction of active electrons decreases with increasing rj^ and for r/g > r/o ~ 4/3 the 
electric field stabilizes the electrons, in that the fraction of active electrons be- 
comes smaller than the one corresponding to equilibria with E = 0. In addition, 
E has similar stabilizing effects on electron NEPs in stellaratorlike equilibria with 
pressure profiles identical to those of tokamaklikc-cquilibria, while it results in an 
increase of the fraction of active ions in rcvcrscd-ficld-pinchlikc equilibria. The 
present results indicate that the radial electric field reduces the NEPs activity 
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in the edge region of tokamaks and stellarators, the reduction of electron NEPs 
being more pronounced than that of ion NEPs. 
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I. INTRODUCTION 



Negative-energy waves are potentially dangerous because they can lead to either 
linear instability |I| or nonlinear, explosive instability - Expressions for 
the second variation of the free energy were derived by Pfirsch and Morrison 
[0] for arbitrary perturbations of general equilibria within the framework of dissi- 
pationless Maxwell- Vlasov and drift kinetic theories. It was also also found that 
negative-energy perturbations exist in any Maxwell- Vlasov equilibrium whenever 
the unperturbed distribution function /^°) of any particle species u deviates from 
monotonicity and/or isotropy in the vicinity of a single point, i. e., whenever the 
condition 

holds (in the frame of reference of minimum equilibrium energy) for any particle 
species u for some position vector x and velocity v and for some local vector k. 
The proof of this result was based on infinitely strongly localized perturbations, 
which correspond to |k| oo. This raises the question of the degree of local- 
ization actually required for NEPs to exist in a certain equilibrium. Studying 
a homogeneous Maxwell- Vlasov plasma P], force-free equilibria with a sheared 
magnetic filed |^ and general one- and two-dimensional equilibria of magnetically 
confined plasmas Correa-Restrepo and Pfirsch showed that NEPs exist 

for any deviation of the equilibrium distribution function of any of the species 
from monotonicity and/or isotropy, without having to impose any restricting 
conditions on k. 

NEPs which are not strongly localized can be investigated more conveniently 
in the framework of Maxwell-drift kinetic theory which eliminates from the outset 
all perturbations with perpendicular wavelengths smaller than the gyro-radius. 
In the context of this theory, for a homogeneous magnetized plasma it was found 
that NEPs exist for any wave vector k with a non-vanishing component parallel 
to the magnetic field (parallel and oblique modes) whenever the condition 

is satisfied for the equilibrium guiding center distribution function /^°) for some 
particle species u and parallel velocity v\\ in the frame of lowest equilibrium en- 
ergy 0]. For the more interesting cases of inhomogeneous magnetically confined 
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plasmas and equilibria depending on just one Cartesian coordinate x and 
cylindrical equilibria with vanishing electric fields [l^, 1^ in addition to parallel 



and oblique modes for which condition (1) also applies, perpendicular NEPs are 
possible. The latter are the most important perturbations because they can ex- 
ist even if v\\{df'^^/dv\\) < 0, which is satisfied e.g. for Maxwellian distribution 
functions for all fii. In plane geometry the pertinent condition is 



< 0, 



dx dx 

where P'-") is the equilibrium plasma pressure. For tokamaklike equilibria with 
singly peaked pressure profiles the existence of both ion and electron perpen- 
dicular NEPs is associated with the critical value 2/3 of the quantity rjy = 
dlnT^/dliaN^ {T^, is the temperature and the density of particle species v) 
which usually governs the onset of the temperature gradient driven modes. For 
cylindrical equilibria an additional regime of NEPs exists, related to the curvature 
of the poloidal magnetic field. Also, for the case of cold-ion equilibria (Tj = 0) a 
large portion of electrons is associated with NEPs (active particles). 

The purpose of the present paper is two-fold: (a) to investigate the impact of 
a radial electric field on NEPs in cylindrical equilibria of magnetically confined 
plasmas and (b) to extend the study to equilibria with Tj ^ 0. The method of 
investigation consists in evaluating the general expression F^"^^ for the second- 
order perturbation energy within the framework of the linearized dissipationless 
Maxwell-drift kinetic theory. This is the subject of Sec. II. The conditions for 
the existence of NEPs are obtained in Sec. III. It turns out that for parallel 
and oblique perturbations condition remains valid, while the condition for 
perpendicular NEPs (which remain the most important perturbations) is mod- 
ified. To apply the condition for perpendicular NEPs in equilibria of magnetic 
confinement systems the equilibrium equations are needed, which are derived in 
Sec. IV. Shearless stellaratorlike equilibria are possible with local Maxwellian 
distribution functions, while tokamaklike and reversed-field-pinchlike equilibria 
can be obtained from shifted Maxwellian distribution functions, which imply net 
toroidal currents. For these kinds of distribution functions and H mode parame- 
ters pertaining to the plasma edge the condition for the existence of perpendicular 
NEPs is applied in Sec. V, and the effect of E on the threshold value of 7],^ is 
examined. In Sec. VI the impact of E on the fraction of active particles is in- 
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vestigated for shear less stellaratorlike, tokamaklike and reversed-field pinchlike 
analytic equilibria. Our conclusions are summarized in Sec. VII. 



II. EQUILIBRIUM AND SECOND-ORDER PERTURBATION 

ENERGY 

We start with a brief outline of the Maxwell-drift kinetic theory adapted to 
the needs of the present study. More details can be found in Ref. [0. 

The expression for the free energy F^"^^ upon arbitrary perturbations of general 
equilibria is given by 

F(2) = / d'x T!i'^\ (3) 



where Tq^^^ is the energy component of the second-order energy-momentum ten- 
sor Tj^^^. To derive the tensor Tj^^^, Pfirsch and Morrison |0] used the following 
modified Hamilton- Jacobi formalism. Let Hu{pi, g', t) be the Hamiltonian for par- 
ticles of species z/ for the perturbed state in a phase space pi, ...,p4, q^, ...,q^ where 
i.l^i'f'i'f) are generalized coordinates so that x = :s.{q^ , q"^ , q^) and correspond- 
ingly p = p{pi,P2,P3), where x is the position vector in normal space; P4,q'^ is 
an additional pair of canonical variables which is needed to describe guiding cen- 
ter motion. Let //^^^(Pj, Q*) be the equilibrium Hamiltonian in the phase space 
Pi, . . . , P4, Q^, ...Q^, and let Su{Pi, q\ t) be a mixed-variable generating function 
for a canonical transformation between pj, and Pi,Q^. The x, t dependence 
of Hi, is given via electromagnetic potentials 0(x, t) and A(x, t), the electric and 
magnetic fields E(x, t) and B(x, t) and their derivatives. The quantities pi and 
are obtained from Si, as 

and 5*1, must be the solution of the modified Hamilton- Jacobi equation 

The time- independent, zeroth-order solution S^^ of Eq. (|), needed to obtain 
Tj;"^^^, is then simply given by the identity transformation S^^ = Y^uPiff- 
The theory can be derived from the Lagrangian 



L = dqdPip,{P,,q\t) 
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-— / d^x(E' 
Stt 



B2). 



(6) 



Here, dqdP = dq^ ■ ■ ■ dq"^ dPi ■ ■ ■ dP^, (p^ are density functions related to the par- 
ticle distribution functions fy and the latter are related to the guiding center 
distribution functions fgi, by Eq. ([T^ ) below. The energy momentum tensor can 
be obtained by using the Euler-Lagrange equations resulting from the variational 
principle 

/.to 

(7) 



5 I Ldt = 0, 



(with {p^, Sy, (j) and A the quantities to be varied) and Noether's theorem. In 
the context of the linearized theory one obtains 



P 



J2 / dqdP 



dqp 



(0) 



\ dq'^ 



dPxdP. 



if(0)p(l) 
' J u Tcr 



dPxdFra 



-^PS E / dqdP 



(0) 
u 



dPi ) , dPx 
dq''^ c '' 



dP.dFlS 



+5i (e / dqdPfl'\n^^^ - ^r^^)) + ^F^'^F^''-') ■ (8) 

Here, the superscripts (0), (1) and (2) denote equilibrium, first- and second-order 
quantities; the tilde signifies that the time is included, i.e. 



(g-) = (g",...,g^) = (ct,x,g^), 

{Pi) = iPo, ■■■,Pa) = {Po, P,Pa), cpo - 
m = {Q',...,Q') = {ct,^,Q') 



dt 



(Pi) = (Po, • • • , Pi) = (Po, P, Pi), Po = constant, 
T-Lyipi.q') = cpo + H^{pi, . . . ,p4,g\ . . . 

Q') = cPo + H,{Pi, ...,P,,Q\...,Q% 
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dqdP = dq^ dPi ■ ■ ■ dP^; = (— 0, A) with A4 = 0; F^i, is the electromag- 
netic tensor; the symbol C*^ signifies covariant derivative of the vector C with 
contravariant components C*: 

where F*, are the Christoffel symbols; the scalar quantity fl^^ — ^ , which 

V dPi y ^ 

replaces ( j in Eq. (46) of Ref. 0], results from the contraction in 



the tensor ( ) . 

V dPj ^ 



The Hamiltonians H^^ for the guiding center motion of a particle species i/, 
which appear in F^'^\ are obtained from Littlejohn's Lagrangian formulation of 



the guiding center theory in the form given by Wimmel ||2^ : 

c 



a: • X - e,<Pl (9) 



with 



= A + g^b, 



ve = c- 



2 

E X B 
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1 (9A B 
E = -V0--— , B = VxA, b = -. 

c at -B 

The Euler-Lagrange equations yield q'^ = \ - h = and the guiding center velocity 
X = V = Vg and as functions of t, x and g^: 

i = V = V,. (t, X, g^) = I^B: + -^E: X b (10) 

and 



g^ = V^(t,x,g^)=^-^E:-B/. (11) 

Here, E* = -V0* - h^^^ B^* = V x A* and = B;^ ■ b. The momenta 
canonically conjugated to x and g^ follow from (^: 

P = = P4 = -TTTT = 0. 12 

ax c ag* 
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Since Eqs. (|T2|) do not contain x and q^, they are constraints between the mo- 
menta and the coordinates. It therefore follows that Hamilton's equations based 
on the usual Hamiltonians corresponding to the above non-standard Lagrangians 
are not the equations of motion. To overcome this difficulty, Dirac's theory of 
constrained dynamics pSf is applied, which yields the Dirac Hamiltonians 



Hu = e,^t + vg, ■ (p - {e,/c)Al) + VW (13) 

Particular solutions of the equations of motion following from the Hamiltonians 
(|13D are the constraints (|T^) . The distribution functions /,^(x, g^, p,p4, t) must 
guarantee that these constraints are satisfied. As concerns this requirement, it 
is important to note that p — A* = and p4 = do not represent special 
values of some constants of motion. Therefore, 5-functions of the constraints are 
not constants of motion either. But f,^ must be proportional to such 5-functions 
and, at the same time, also a constant of motion. Both conditions are uniquely 
satisfied by 

/. = Sip,)6 (p - ^Al^ Bt^^f,, (x, , (14) 

where the guiding center distribution functions fg^ are constants of motion and 
solutions of the drift kinetic differential equations 

In the present paper cylindrical equilibria are considered. With the coordi- 
nates q^,q^,q^ specified to be the cylindrical coordinates r, 0, z with unit basis 
vectors e^, e^, the equilibrium vector potential and magnetic field are given by 

= 4°^(r)e9 + 4°)(r)e, (16) 

and 

BW = i?(°)(r)e, + i?f(r)e„ (17) 

with 

\rA^')' = Bf\ {Af^)' = -Bf\ (18) 
r 

where the prime denotes differentiation with respect to r. The equilibrium electric 
field can be expressed in terms of the scalar potential 0*^°)(r) as 

= -V0(°) = -(0(°))'e,. (19) 
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For the equilibria defined above, the guiding center velocity [Eq. (]T0|)] becomes 



|b(0) + v 



(0) 



*(o) 



+ 



m 



e, X b(o)) , 
(20) 



where v^'' = b*^°) x {yfj x b'^'^^^ is the perpendicular component of v^^^ consisting 



of the E X B, V-B, curvature, and polarization drifts; 



*{0) 



B 



*(o) 



b(o) X (B;^(°) X b(°)). For thermal particles it holds that Blr ^ and 



B 



*(0) 



B^(o) ■ b(o) and 

^(0) 



~ O{rgi,/ro) where r^,^ is the thermal Larmor radius for the particle species 

V and To the macroscopic scale length. v^°^ has no r-component and therefore r is 
a constant of motion. Since there is also no force parallel to B'^°\ another constant 
of motion is the parallel guiding center velocity f y . The guiding center distribution 
functions f'f^ are therefore functions of r, v\\ and the magnetic moment /i. From 
Eq. (^) it follows that Vf^ = and hence the Dirac Hamiltonians [Eq. (|T^] 
are written in the form 



(21) 



The general expression for the second order perturbation energy [Eq.(^] is 
evaluated for these equilibria and for initial perturbations A^^) = and A^^^ = 0. 
It is also shown a posteriori that one can choose initial perturbations such that 
the charge density p*^^^ vanishes without changing the particle contributions to 
the energy. Thus, choosing perturbations of this kind, we can set from the outset 



7(1) 



AW 







After a lengthy derivation, which can be conducted along the lines of that for 
cylindrical equilibria with E = reported in detail in Appendix A of Ref. |19 
is cast in the concise form: 



ir(2) 



/ S{r)drdv\\d^ 
^ ^ dv\\ 



B 



*{0) 



k- V 



dr 



(22) 
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Here, u^'^^^ = {e^Bl^^ )/{c'my), G^J\r,q^, fi) are arbitrary first-order quantities 
relating to the generating functions S*^^-* for the perturbations; h*^^^ = B**^'^) / B*^^^ 
; k = keee + kze^ is the wave vector lying in magnetic surfaces; k\\ and k± are its 
components parallel and perpendicular to B*^"), respectively, and 



2tt , I 27r 



S{r) = r / « / ^ dOdz 

is a normalization surface, where 6q and zq are constants. We note that F*^^) 
depends on G\^^ only via IG^^^^p. The first-order charge density p*^^-* is a 
integral over an expression that is linear in G^J^^ . One can therefore satisfy the 
relation p^^^ = by a proper distribution of positive and negative values of G^,^-* 
on which F^'^^ does not depend. 

Compared with the corresponding expression for equilibria with E'^'') = 
(Eq. (37) of Ref. [|l^), -F*-^-* contains new terms stemming from b*(°^ and from 
the E X B and polarization drift components of v^j^. In particular, as shown in 
Section III the E x B drift modifies the condition for the existence of NEPs with 
wave vectors perpendicular to B*-*^-*. 

III. CONDITIONS FOR THE EXISTENCE OF 
NEGATIVE-ENERGY PERTURBATIONS 

The derivations in this section are very similar to those concerning equilibria 



with E = |T^, [T^, |T9[ so that details need not be given here. 

The following conditions must be satisfied only locally in r, v\\ and /x, and 
refer to the frame of reference of minimum energy. 

Parallel perturbations (k±_ = 0) NEPs exist when 

> (23) 



dv\\ 



is satisfied for at least one particle species. This condition, which was first de- 
rived by Pfirsch and Morrison for a homogeneous magnetized plasma [^, agrees 
with those obtained by Correa-Restrepo and Pfirsch for several Maxwell- Vlasov 
equilibria 

Oblique perturbations {k\\ ^ and k_[_ ^ 0) If condition ( p3D is satisfied for at 
least one particle species only perturbations with wave vectors satisfying in 
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addition the relations 



with 



and 



fell fell 

< min(A^, M^) or ^ > max(A^, M^), (24) 

K I Ki 



(0) 

7)11 V 



1 



^.(0) 5^ ^wo) 5^11 j y 5^11 y 

can have negative energy. The orders of magnitude of A^, and Mj^ depend on the 
particle energy. For example, if 

m^vl ^ fiB^^^ ^ \e^<p\ ^ (25) 

with (f){oo) = it holds that 

lAJ ^ IMJ ^ ^ « 1. (26) 

Relation (|26|) indicates that condition (|2^) imposes no essential restriction on the 
magnitude or the orientation of the wave vectors associated with NEPs. 
If 

-11^ <0' (27) 

a condition which is satisfied at all points of a Maxwellian distribution function, 
NEPs exist if, in addition to (ETD, it holds that 



fell 

min(A,„M^) < < max(A,„ M^). (28) 

K± 



For the scaling (P3|) condition (53) implies that 



^ ~ !> « 1, (29) 
k± ro 

which indicates that the most important NEPs, in the sense that the less restric- 
tive condition ( jn\ ) is involved, are associated with nearly perpendicular wave 
vectors. It may be noted that for a homogeneous magnetized plasma in thermal 
equilibrium, although condition (^) is satisfied, NEPs are not possible because 
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My = Aj, = and therefore condition (^) is not satisfied. This also follows 
from condition (^) which is pertinent for the existence of NEPs in a homoge- 
neous magnetized plasma and is not satisfied in thermal equilibrium, i.e., for 
Maxwellian distribution functions. 



Perpendicular perturbations {k\\ = 0) In this case the second-order perturbation 
energy [Eq. (p2D] reduces to 



?*(0) Tjr / , \ 2 



with 



and 



(B^^^y \eu ^(0) iu{o)\i ( (0) 2V' 



(31) 



- 9r 5^(0) ' ^^^^ 

where Wj,± = fxB^^^ is the perpendicular particle energy. Eq. (|30D implies that 
< for any k± whenever the condition 

Z,Q, < (33) 

is satisfied irrespective of the sign of v\\{df'^^ / dv\\). Therefore, there are two 
regimes of NEPs which are determined by the relations 

< and > (34) 

and 

> and Qy < 0. (35) 

For the evaluation of conditions ([3^ ) and (^) the equilibrium equations are re- 
quired, which are constructed in the following section. 

IV. QUASINEUTRAL MAXWELL-DRIFT KINETIC 
EQUILIBRIUM EQUATIONS 

The equilibria must satisfy 

V ■ E(°) = 47rp(°) (36) 
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and 



c 



(37) 



where the charge density p*-*^-* and current density j*-"-* are expressed self-consistently 
in terms of the guiding center distributions functions f^'^^ in the context of the 
Maxwell-drift kinetic theory (see Eqs. (8.14) and (8.15) of Ref. [^). For the 
system under consideration, owing to the presence of 'E^^\ the set of equilibrium 
equations following from Eqs. (3B) and ( P7| ) are rather complicated. For this rea- 
son we employ the quasineutral Maxwell-drift kinetic theory which can be derived 
self-consistently by dropping the electric- field-energy term in the Lagrangian (|]). 
A similar method was employed in Refs. 0, and [^. Consequently, Eq. 
(|36D is replaced by the quasineutrality condition, which is explicitly given by 



(0) .(0)^/(0) _ 



91^ ^^fi^ ^ 



(0)^ 



:b(°) =0 (38) 



and Eq. (|37|) by 



curl / dv\\djj,B*^^^ 

1, 



(0) f (0) 
9^ 



V 



(E(o))2 



,(0) 
' gv 



vi?)) X + 



2m 



fgv^ 



(0) 
91^ 



—V X B(°) 

4:71 



(39) 



The first terms on the left-hand sides of Eqs. (^) and (^) represent guiding- 
center charge and current density contributions, respectively, and the other terms 
polarization and magnetization contributions. 
We consider equilibria of the following kind: 



1. The distribution functions are specified to be local shifted Maxwellians 



(0) 



271 



fiB^^\r) + l/2m^ 



Ti°)(r) 3/2 



exp 



v;(o)(r 



Ti°)(r) 



(40) 

where Nj)^^ and Tj°^ are, respectively, the density and temperature for par- 
ticles of species u. They can describe cylindrical tokamaklike, reversed-field 
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pinchlike and for Vj''-' = shearless stellaratorlike plasmas, which are close 
to thermal equilibrium. For the former equilibria the shift velocities Vj''-' 
satisfy 

- ^ < 1 (41) 



{vy)th ro 

and, as shown later, lead to a non- vanishing toroidal current. 

2. Since a radial electric field may play a role in the L-H transition of magnetic 
confinement systems, e.g. [^, |2^, we adopt for the ion electrostatic energy 
the scaling 



.(0) 
i 

which is satisfied in the edge region. 



6^0-' I ^Tf^ 0(°)(oo) = O, (42) 



Using the above assumptions, neglecting small terms of the order Vg^/rQ and 
suppressing the superscript (0) from the equilibrium quantities, Eq. (|38D and the 
9- and z- components of Eq. (|39|), respectively, yield 

J2e,N, = 0, (43) 

V 

jg = bge,N,{Vi - Ve) + c^P' = -^B[ (44) 

B An 



and 



where 



Jz = KeM{V, - Ve) - c-^P' = -^-{rBeY, (45) 

B 4:71 r 



P = Y.I dv\^dfiBt\\fiBfg, = (46) 
For Vu = for all u Eqs. (|4^) and (^5|), respectively, reduce to 

hp' = (47) 

B An ^ ' 

and 

The solutions of Eqs. (|47|) and ( ^Sf ) satisfy the relation .Bg = a-p, with a = 
constant. They are singular at r = and therefore a = 0. For Bg = Eq. (^) is 
satisfied identically and the only possible equilibrium, which is described by Eq. 
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(^7[), is a 6'-pinch or shearless stellaratorlike configuration with vanishing axial 
current, which is examined in Sec. VI. 




Multiplying Eqs. (44) and ( ^5|) by the integrating factors Bz and Bg, re- 
spectively, and adding the resulting equations, one obtains the pressure balance 
relation 

|:(^ + ^ 1+^ = 0, (49) 

which will be used in place of Eq. 

Summarizing, quasineutral equilibria can be described by the set of Eqs. (^31), 
(PI), (^) and (|49|). Four out of the eight functions involved must be assigned, 
e.g., P(r), Bz{r), the shift velocity difference Vi{r) — Ve{r) and Tj(r); then, Bg{r) 
can be obtained from Eq. (^), Ni{r) from Eq. (0), Ne{r) from Eq. ( ^3| ) and 
Te(r) from Eq. (|^). Analytic solutions, which are required for determining 
the portion of active particles in equilibria of magnetic confinement systems, are 
constructed in Sec. VI. 

V. PERPENDICULAR NEPs IN EQUILIBRIA OF MAGNETIC 
CONFINEMENT SYSTEMS 

In this section condition ( |33D for the existence of perpendicular NEPs is applied 
to the equilibria defined in . II and IV. For distribution functions of the form 
( plD and the pressure balance relation (^), the quantity Qu [Eq. (^)] reduces 



to 

Qu = jfU, (50) 

where 

2 

Uu = I - - + Vu^lu + e2u, (51) 

e2u = (53) 

and 

Depending on the value of Ti/T^, the effect of E on perpendicular NEPs is exam- 
ined in the following two regions. 
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Assuming the scaling (^) to hold it can be shown that the pressure gradient 
and a term relating to the curvature of Be dominate in Zi^, i.e. 

1 + 2^] =K. (55) 




W, 



Condition ^ 

Rv^-ttUu < 0. (56) 



Relation ( p6D is identical to the corresponding one in equilibria with E = (Re- 
lation (58) of Ref. [0). For singly peaked density and the temperature profiles 
and therefore rjy > {] for all which is the most common case in equilibria of 
magnetic confinement systems, there are two regimes of NEPs depending on the 
sign of Ry [Eq. 



a) -R^ < Condition (|56D implies that Uy <Q must hold. The last two terms of 
Uy [Eqs. (0)-(0)] become non-negative and vanish for Wy\\ and Wy± 0. 
Consequently, Uu < is satisfied whenever 

Vu > I (57) 
The existence of perpendicular ion NEPs for any kj_ is therefore related to the 



threshold value of 2/3 of the quantity 7]^. As discussed in Ref. this threshold 
value is sub-critical in the sense that it is lower than the critical value rj'l 1 for 
linear stability of temperature-gradient-driven modes. 

b) Ry > Q The condition for the existence of perpendicular NEPs becomes 
Uy > 0. In this case the quantity rj^eiy + can be either positive or negative 
and therefore no restriction is imposed on rjy. It may be noted that for plane 
equilibria it holds that Ry = P' < and therefore the second regime of NEPs is 
associated with the curvature of Be- 

B) |f>/?e 
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If the scaling (|i2| ) holds the term Cycf)' related to the E x B drift dominates 
in [Eq. (0)], i.e. 

Condition fB^) becomes then 



e^'-^U, > 0. (59) 

Relation ( |59D shows that the existence of perpendicular NEPs depends on the 
sign of the particle species charge and the polarity of E. Henceforth and up to 
the subsection "Reversed-field-pinchlike equilibria" E 7^ will refer to this case 
(Tj/Tg > Pc)- In the edge region the radial electric field is usually negative |2^, ^ 



It is noted here that the impact of the polarity of an externally induced radial 
electric field was investigated experimentally It was found that whereas the 
energy confinement in H-modes with E > is at least as good as in those with 
E < 0, the ratio of the ion confinement time to the energy confinement time is 
about three times lower in the former case. We examine therefore in the following 
NEPs for ions and electrons in equilibria with 0' > 0. 



Ions In this case eu(f)' is positive and condition (^) is satisfied whenever Ui < 0. 
Depending on the sign of Ri, there are two regimes of NEPs: a) If -Rj < 0, t/j < 
is satisfied whenever rji > 2/3 and b) if Ri > no restriction is imposed on rji. It 
is pointed out, however, that for E = the condition associated with this second 
regime is t/j > 0. As shown in Sec. VI, this difference affects the fraction of active 
ions. 



Electrons Condition (59) is satisfied whenever f/e > 0. This yields 

2 , 

r^e < - (1 + Ve^le + (^2e) ■ 



(60) 



For cold electrons (We|| ~^ and IVe± ~^ 0) condition ( pDj ) implies that NEPs exist 
whenever ?7e < 2/3. This indicates that E has a "stabilizing" effect on electron 
NEPs for large values of ?7e. Owing to hot electrons however, condition (60) does 
not yield an upper threshold value of ?7e because electrons with non-vanishing 
energy activate NEPs in the regime where > 2/3. To determine the value of 
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rje for which half of the electrons are active, condition ( |60D is written in the form 

C.„.)i^ + Z,.M^^>|-i (61) 



where 



2 52 T, 



^e(r)= 1 + --^^ and D,(r) 



47r / „, B, 



2 



1 + ^^ P' + 



52 Ti \ 47rr 



(62) 



For a magnetic confinement system it holds that Bg/B"^ ^ P/{B^/87t) = (3 with 
max/5 0.1 and therefore Ce ~ -De ~ 1- Consequently, condition (^) implies 
that nearly half of the electrons are active whenever it holds that 3/2 — l/?7e ~ 3/4. 
This yields 

t ^ 4/3. (63) 

Therefore, if E 7^ less than half of the electrons are active whenever rjf. > 77° 
and this portion decreases as rje takes larger values. On the other side, if E = 
more than half of the electrons are active whenever rj^ > 77° with this portion 
increasing as rje takes larger values. 

VI. ANALYTIC EQUILIBRIUM SOLUTIONS AND ACTIVE 

PARTICLES 

In this section the portion of active particles is determined on the basis of ana- 
lytic shearless stellar atorlike, tokamaklike and reversed field pinchlike equilibrium 
solutions. 

Shearless stellaratorlike (6 -pinch) equilibria 
We consider the following profiles: 

P = P(0)(l-p2), (64) 

Ni = A^i(0)(l — and Ti = Tj(0)(l — p^Y'^, where p = r /tq and tq correspond- 
ing to the plasma surface. Eqs. ( ^9]) with Bq = 0, J2u ^vNy = and P = J2u ^uTu 
yield then 



B.. = B 



zO 



l-/3o(l-p') (65) 



] = yBz X = po ttt e^, 

^o'^°[l-/?o(l-p^)]'/' 
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= A^e(0)(l -p2)« and = re(0)(l -p^)!-?. Here, B,o is the external constant 
"toroidal" magnetic field , /3o = P{Q)/{B'^/M) and the parameter ^ (0 < ^ < 1) 
determines equilibria with different values of r]^, i.e., 

_dwr^ i-e .... 
= 5h^ = 

Ion and electron NEPs are now examined separately. 



Ions Since the magnetic field lines are straight, it holds that Ri = P' < and 
therefore ion NEPs exist only in equilibria with rji > 2/3. The pertinent condition 
Ui < becomes 

Ti [ 2(1 - ^) J Ti 2 ?7i 

Relation ( |67D implies that: 

1. The portion of active ions increases as rii takes larger values. In particular, 
for a fiat temperature and peaked density profile there are no active ions; 
for rji = 1 one third of the thermal ions are active and for r]i = 2 this fraction 
becomes 2/3; for a flat density and a peaked temperature profile {rji oo) 
all ions are active; 

2. The portion of active ions increases from the center p = to the edge region 
p = l. 

Electrons For E = the situation is similar to the foregoing one for ions. For 
E 7^ the condition for the existence of electron NEPs is ?7e > and therefore 
the fractions of active electrons and ions are complementary to each other. Thus, 
as also discussed in Sec. V, the electric field stabilizes electron NEPs for rje > 4/3, 
e.g., for the equilibrium profiles (|64D -(|66D 1/3 of the thermal electrons are active 
when r/e = 2, while the corresponding fraction for the equilibrium with E = is 
2/3. In addition, the fraction of active electrons decreases from the center to the 
edge. This indicates that in the presence of E self-sustained turbulence associated 
with electron NEPs should be reduced in the edge. 
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Tokamaklike (screw pinch) equilibria 
The following profiles correspond to a special solution of Eq. (|49| ) : 



S^(0) + 87rP(0)(l -a2)p^ 



1/2 



(68) 



a is a parameter which describes the magnetic properties of the plasma, i.e. the 
plasma is diamagnetic for a < 1 and paramagnetic for a > 1; 



Be = 2J7rP{0)ap; 



(69) 



constant axial current density; N^, = A^,^(0)(1 — p^)^ and T^, = Tj,(0)(l — p'^Y ^ 
with u = i,e. 

Ion and electrons NEPs are now examined for t]^ = 1, which is close to linear 
stability threshold for gradient temperature driven modes. 



Ions For E = the portion of active ions is determined by conditions (|^) and 
(PSD which respectively become 



Will 1/2 ^ , 

< - — - 1 and 



and 

Wiw 1/2 
— — -1) and 
Wi± 2 



l-\(3{l-p'){a'-2) 



l-pa\l-p^ 



T,. 



l-l3a\l-p^ 



Wi\\ 1 

~T~ ^2 
(70) 



T,; 2 



(71) 



Relations (^) and (|7TD imply that: 



1. The portion of active ions increases with a, i.e., it is smaller in a paramag- 
netic and larger in a diamagnetic system. The particular cases of a strongly 
diamagnetic plasma (a — > 0), of an equilibrium with constant "toroidal" 
magnetic field (a^ = 1) and of a paramagnetic plasma {a^ = 2) are illus- 
trated, in Figs. 1, 2 and 3, respectively. The fractions of active ions are 
nearly 1/3 for a ^ 0, 1/2 for = 1 and 2/3 for = 2. It is noted that for 
a — i> only the branch ([70|), associated with the threshold value rji = 2/3 
contributes while for = 2 exclusively the branch (|71|) associated with the 
curvature of the poloidal field hues contributes. 
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FIG. 1. The portion of active ions for a strongly diamagnetic plasma with 
E = and rji = 1 which is deduced from Eq. ([70| ) [ao(p) = 1 + — p^)]- 
The dotted area stands for the active particles at the center (p = 0), while 
the area filled by circles for the additional active particles at the edge (p = 

1)- 



2. In all regimes the fraction of active ions increases from the center to the 
edge. In Figs. 1, 2 and 3 the dotted area stands for the active particles at 
the center (p = 0), while the area filled by circles for the additional active 
particles at the edge (p = 1). 

It is noted here that for E = similar results hold for electrons. 
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FIG. 2. The portion of active ions for the equihbrium with E = 0, r/j = 1 and 
Bz = constant, which is deduced from Eqs. (ffOD and (|7lD [ai(p) = 1 + (/9/2)(l — 
p% 6i(p) = 1-/3(1- p2)]. 




FIG. 3. The portion of active ions for the equilibrium of a paramagnetic plasma 
with E = and r/e = 1, which is deduced from Eq. (|7l]) [b2{p) = 1 - 2/5(1 - p^)]. 
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For E 7^ 0, active ions obtain from condition Ui < (irrespective of the sign 
of Ri) which leads to 



+ 



Wi, 1 



(72) 




FIG. 4. The portion of active ions for the equihbrium with E 7^ 0, ?7j = 1 and Bz = 
constant, which is deduced from Eq. (|72D [a^^p) = 1 + (/9/2)(l — p^), &4(p) = 1 — 
/3(1 — p^)]. The excess portion at the edge indicated by circles nearly compensates 
for the excess portion at the center indicated by stars. 



Relation ( [72D implies that: 

1. The portion of active ions is nearly independent of the magnetic properties 
of the plasma; it is approximately 1/3 for any value of a. 

2. The portion of active ions can either be nearly independent of p, e.g., for 
an equilibrium with constant Bz (a^ = 1) (Fig. 4) or decreases from the 
center to the edge, e.g. for a paramagnetic plasma = 2 (Fig. 5), while 
this portion always increases for equilibria with E = 0. 



Thus, E leads to a reduction of active ions. 
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FIG. 5. The portion of active ions for a strongly diamagnetic plasma with E 7^ 
and rji = 1 deduced from Eq. ( [72D [05 (p) = 1 — 2/5(1 — p^)]. The area filled by 
stars represents the excess portion at the plasma center. 



Electrons Recalling that the portion of active electrons is the same as that of 
active ions when E = 0, and complementary when E 7^ 0, respectively, the former 
portion can be determined on the basis of the foregoing analysis for ions. Thus, in 
addition to the stabilizing effect of E for r/g > 4/3, the fraction of active electrons 
(a) becomes nearly independent of the magnetic properties of the plasma and (b) 
can decrease from the center to the edge, e.g., for the most common case of a 
diamagnetic plasma. 

Reversed-field-pinchlike (force-free) equilibria 

The solution of Eq. with P' = leads to B, = B,{0)Jo{p) and Eg = 

Bz{0)Ji{p), where Jq and Ji are Bessel functions. These profiles satisfactorily 
describe the central region of the relaxed state of a reversed-field pinch By 
appropriately assigning Vi{r) — V^(r), one can derive equilibria with a variety of 
density and temperature profiles for which NEPs exist and a considerable fraction 
of active ions and electrons are involved. From the equilibria considered it turns 
out that E (a) does not affect the electron NEPs and (b) enhances the fraction 
of active ions. 
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As an example, we consider the most common equilibrium with constant den- 
sity and temperature profiles: 



V 



V 



T,. 



(73) 



For E = 0, with the aid of relation ( |55D condition ( |33D becomes 




(74) 



for any particle species u. Therefore there are neither ion nor electron NEPs. 
If E 7^ 0, NEPs exist whenever the condition 



particle species charge in condition ([75l), for > all ions are active, while the 
active electrons are not affected. 



The impact of a radial electric field on negative-energy perturbations (NEPs) in 
cylindrical equilibria of magnetically confined plasmas was investigated within 
the framework of linearized dissipationless Maxwell-drift kinetic theory. The in- 
vestigation consisted in evaluating the general expression for the second-order 
perturbation energy derived by Pfirsch and Morrison for the equilibria under 
consideration and for vanishing initial field perturbations; then the conditions for 
the existence of NEPs were obtained. 

The electric field E does not affect the following condition for perturbations 
with wave vectors parallel and oblique to the equilibrium magnetic field (fcy ^ 0): 
If the equilibrium guiding center distribution function (r, v\\ , /i) of at least one 
particle species u satisfies the relation U||((9/^°V^'^||) > locally in r,v\\ and /i, 
parallel and oblique NEPs exist with no essential restriction on k. The condition 
for the existence of perpendicular NEPs {k\\ = 0), which holds regardless of the 
sign of v\\{dfj^^ydv\\), is modified. For |ej0| ~ Tj the effect of E on perpendicular 
NEPs depends on the value of T^/Te, i.e., a) for T,/Te < /3c ~ P/{B^/87t) the 
electric field has no effect, and b) for Tj/Tg > /5c, a case which is of operational 




(75) 



following from relations ( pSj) and (|58|), is satisfied. Owing to the presence of the 



VII. CONCLUSIONS 
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interest in magnetic confinement systems, the existence of perpendicular NEPs 
depends on tlie sign of tlie particle species charge and the polarity of E [Relation 
(^]. For E < it was found that: 



1. For cylindrical tokamaklike equilibria described by local shifted Maxwellian 

distribution functions and singly peaked pressure profiles there exist two 

regimes of NEPs for both ions and electrons. One regime is associated 

with the curvature of the poloidal magnetic field. In the other regime the 

d In T 

threshold value 2/3 of r^j = ^i^jy is involved for ion NEPs, as in equilibria 
with E = 0, while a critical value of 7]^ does not occur for the existence of 
electron NEPs. However, E has the following "stabilizing" effects on both 
particle species: 

• The portion of particles associated with NEPs (active particles) is 
nearly independent of the plasma magnetic properties, i.e. it is nearly 
the same in a diamagnetic and in a paramagnetic plasma, while in 
equilibria with E = this portion is much larger in a paramagnetic 
than in a diamagnetic plasma. 

• The portion of active particles can be either constant or decreases 
from the center to the edge, e.g., in the case of active electrons of 
a diamagnetic plasma, while it always increases in the corresponding 
equilibria with E = 0. 

In particular, the fraction of active electrons decreases with increasing rje 
and for rje > t]o ~ 4/3 the electric field stabilizes electron NEPs in the 
sense that the fraction of active electrons becomes smaller than the one 
corresponding to equilibria with E = 0. 

2. In shearless stellaratorlike equilibria described by local Maxwellian distri- 
bution functions and pressure profiles identical to those of tokamaklike- 
equilibria, E leads to similar stabilizing effects on electron NEPs; namely, 
it reduces the fraction of active electrons (a) for r^e > r/o ~ 4/3 and (b) from 
the center to the edge. 

In addition, irrespective of the value of Tj/Tg, E does not affect the electron NEPs 
in reversed-field pinchlike equilibria but "destabilizes" the ion NEPs in the sense 
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that it enhances the portion of active ions. For example, for an equihbrium with 
constant density and temperature profiles all ions are active in the presence of E, 
while there are not active ions when E = 0. 

The present results indicate that a radial electric field leads to a reduction of 
the NEPs activity in the edge region of tokamaks and stellarators. For electrons, 
which may mainly contribute to anomalous transport, this reduction is more 
pronounced. 

Finally, it may be noted that according to the results of our previous p^ , 
pj| and in the present studies, the curvature of the poloidal magnetic field is 
unfavorable in the sense that it gives rise to an increase of NEPs activity. It can be 
speculated that this is true for an arbitrary magnetic field configuration. To check 
this conjecture, it is interesting to investigate NEPs in a toroidal equilibrium, e.g., 
a tokamak, in which the toroidal magnetic field is favorably curved on the inside 
and unfavorably on the outside of the torus. Such a study might also reveal the 
effect of toroidicity on other aspects of NEPs, e.g., the threshold value i],^ = 2/3. 
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